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The  Structure  of  Light  Radiance  in  the  Depth  of  the  Sea 
for  the  Problems  of  Horizontal  Visibility  of  Submerged  Objects 

Vladimir  I.  Haltrin* 
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ABSTRACT 

The  spatial  and  angular  structure  of  light  radiance  in  the  depth  of  seawater  can  be  precisely  computed  with  such 
numerical  programs  like  Hydrolight.  Unfortunately  such  calculations  demand  a  lot  of  computational  time  if  used  for 
problems  of  visibility  of  submerged  objects.  In  this  case,  due  to  the  much  greater  computational  efficiency  and 
logarithmic  dependence  of  visibility  parameters  on  radiance  values,  analytical  radiative  transfer  methods  are  preferable.  In 
this  presentation  a  new  analytical  radiative  transfer  method  is  proposed  to  solve  visibility  problems  in  seawater.  The 
angular  structure  of  depth-dependent  light  radiance  in  the  sea  is  obtained  in  the  framework  of  previously  published  self- 
consistent  approximation.  The  analytical  expressions  proposed  in  this  paper  allow  to  compute  a  horizontal  visibility  of 
submerged  objects.  This  approach  can  be  applied  to  problems  of  diver  visibility  in  shallow  water  bodies  of  arbitrary 
turbidity  and  illumination. 

Keywords:  shallow  water,  radiative  transfer,  ocean  optics,  horizontal  visibility. 


1.  INTRODUCTION 

The  purpose  of  this  paper  is  to  obtain  the  simplest  analytic  solution  to  the  radiative  transfer  problem  in  a  shallow 
water  with  wavy  surface  and  reflecting  bottom  that  takes  into  account  all  orders  of  scattering  inside  water  as  well  as  all 
orders  of  scattering  (or  reflection)  from  bottom  and  sea  surface.  The  approach  used  here  is  a  self-consistent  approach 
(SCA)  [1-3]  that  was  originally  devised  to  obtain  upward  and  downward  irradiances  in  shallow  sea  with  arbitrary 
inherent  optical  properties.  Potentially  SCA  is  capable  to  produce  radiances  of  total  light  distribution  in  shallow  sea 
because  the  main  SCA  radiative  transfer  equation  has  a  source  term  that  depends  only  on  radiance  of  a  source  light  and 
irradiances  of  scattered  light.  This  potentiality  was  utilized  in  this  paper  to  obtain  total  radiances  of  light  inside  shallow 
sea  and  horizontal  irradiances  that  can  be  used  for  estimation  of  horizontal  visibility.  Calculation  of  horizontal  visibility 
is  considered  to  be  the  next  task  and  can  be  easily  performed  using  the  results  obtained  in  this  paper. 


2.  RADIATIVE  TRANSFER  PROBLEM  IN  SHALLOW  WATER  WITH  COMBINED 
ILLUMINATION  BY  SKY  AND  SUN 

We  start  from  the  standard  scalar  radiative  transfer  equation  for  total  radiances  of  light  L,  in  the  ocean  depth: 

+  c)  m&w) m  'fafo*  ^'f-i  ^^(cosy  (0 
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here  c  -  a  +  b  is  an  attenuation  (extinction)  coefficient,  a  is  an  absorption  coefficient,  b  is  a  beam  scattering 
coefficient,  p(cosy)  is  a  scattering  phase  function  of  seawater  normalized  as  0.5 J  ^/?(cosy)siny  dy  —  1,  y  is  a 

scattering  angle,  Z  is  a  Cartesian  depth  coordinate  originated  in  the  water  surface  (z  =  0)  and  going  down  to  the  bottom 
(2  =  %  )  .the  X  -  y  plane  coincides  with  the  water  surface  when  it  is  calm;  the  zenith  angle  6  -  COS  p  is  measured 
from  0 z  axis,  and  the  azimuth  angle  cp  is  measured  from  Ox  axis.  The  scattering  angle  y  is  determined  by  the 
following  formula: 

cosy  =  pp'  +  cos(<p  -cp')  (2) 

here  p  =  COS 6,  pi'  =  COsQ' .  The  source  portion  of  light  just  below  the  sea  surface  consist  of  the  following  two 
components:  1)  diffuse  component  originated  from  the  light  of  the  sky  transformed  by  the  wavy  surface  into  Lambertian 
source,  and  2)  direct  sunlight  penetrated  into  the  sea  water. 

We  consider  the  case  of  shallow  coastal  part  of  the  ocean  with  the  bottom  located  at  depth  zB  and  reflecting  light 

Lambertially  with  bottom  albedo  AB .  The  depth  of  the  ocean  may  be  less  than  one  optical  depth,  1 1 C. 

According  to  previous  investigations  [4]  the  wavy  sea  surface  in  the  first  approximation  may  be  considered  as  a 
Lambertian  reflector.  The  albedo  of  the  surface  A5  for  ascending  light  is  a  function  of  wind  speed,  sea  foam  albedo,  and 

refraction  coefficient  nw  of  seawater. 

We  solve  Eq.  (1)  for  in-water  radiances  using  self-consistent  approach  [1]  to  the  radiative  transfer  of  light  in 
seawater.  This  approach  consists  of  reducing  the  scattering  phase  function  /?(cosy)  to  the  one-parameter  transport  phase 
function,  and  applying  two  experimentally-based  conditions  that  connect  the  average  cosines  of  diffuse  portion  of  light 
in  water: 

/?(cosy)  —*  pT (cosy)  =  2 B  +  2(1  -  2B)<5(1  -  cosy),  6(1  -  cosy)  =  2jv5(cp  -  cp  )d(p  —  p  ),  (3) 

pd=\/(2-p),  ft- 1/(2 +  M).  (4) 

here  B  =  0.5 J*  p(p)  dp  is  a  backscattering  probability,  pd  is  a  downward  average  cosine,  lllt  is  an  upward  average 
cosine,  and  p  is  a  total  average  cosine  [1-3].  The  average  cosines  will  be  defined  through  irradiances  of  light  later.  Here 
and  further  <5(jc)  denotes  a  Dirack’s  delta-function. 

The  self-consistent  approach  has  one  huge  advantage  over  all  other  approximate  theories:  it  is  valid  for  all 
possible  values  of  inherent  optical  properties:  extinction  c,  scattering  b,  and  absorption,  coefficients,  i.e.  it  is  good  for 
0  £  a,b,c  <  oo. 

By  applying  a  self-consistent  condition  to  Eq.  (1)  we  obtain  the  following  equation: 

where  CC  =  a  +  2 bB  is  a  re-normalized  extinction  coefficient,  and  bB  —  bB  is  a  backscattering  coefficient. 

Now,  let  us  split  the  total  light  radiance  Lt  into  two  components:  direct  L ^  and  scattered  L  light, 
L,  =  Lq  +  L .  If  the  penetrated  light  of  the  sky  and  sun  just  below  of  the  sea  surface  is  expressed  as: 

Lq(0,p,cp)  =  L0D+L0S6(p-p0)S(cp-n),  (6) 

where  L0D  is  a  radiance  of  diffuse  light,  Los  is  a  radiance  of  sunlight,  then  the  radiance  of  the  source  will  be: 


Proc.  of  SPIE  Vol.  5544  373 


(7) 


L„{z^<p) 


Lod  e'a  z/,‘  +  Los  e~°  J,"°  6(n-  p0 )d(<p  -  jt), 

0, 


AtS:  0, 

jM<  0, 


here  /i(!  =  cosd0  is  a  cosine  of  the  angle  60  between  direction  of  direct  solar  light  in  the  water  and  axis  0z.  It  is 
determined  through  seawater  refractive  index  nw  and  solar  zenith  angle  z®  via  Snellius  law,  sin  Z©  =  nw  sin  0O,  or 


~t*o~ 


sin  z® 


■<&- 


By  inserting  L,(z, /-*,<?)  =  Lq(z,H,<p)  +  L(Z,fa)  into  Ecl‘  (5)  with  Lq  Siven  by  Ecf  we  have  the  folIowinS 
equation  for  scattered  radiance  of  light  inside  seawater: 


where 


=  bBJ [d[x'  L(z,ti')  +  Q(z), 
-  h  L,  q(a  z)  +  ] ,  *t)  -fan 


is  the  source  function.  . 

To  proceed  further  from  Eqs.  (9)-(10)  let  us  define  diffuse  irradiances  and  average  cosines  as  follows: 

Downward  ( Hd )  and  upward  ( Hu)  scalar  irradiances: 

Hd  (z)  =  d(pJo  dfx  L{z,fa)  =  2jtjQdn  L(z,fa) , 

H„(z)  =  fo  dvf  dfJiL{z,n)  s  2jtf^dfiL(z,fa), 

Downward  ( Ed )  and  upward  ( Eu)  irradiances: 

Ed(z)  =  dcp^ndixLiz,^)  s2jijQndnL(z,^), 

Eu(z)  =  ~J0  dxpj  /J.dfiL(z,^)  =  -2jtJ  fxdfj,L(z,fj) 


(9) 

(10) 


(11) 

(12) 


(13) 

(14) 


Ji<i  ~  Ej/Hj  ,  K-EJHu ,  i u  =  (Ed-Eu)/(Hd+Hd).  (15) 

Using  Eqs.  (1 1)-(15)  we  can  rewrite  Eq.  (9)  as  follows: 

ln—  +  a)  L(z,n)  =  ^H(2-  jd)Ed  (z)  +  (2  +  JI)Eu(z)]  +  Q(z) ,  (16) 

\  dz  I 

where  radiance  of  scattered  light  is  defined  through  the  irradiances  of  that  light  and  source  function  of  external  light 

given  by  Eq.  (10).  . 

Equation  (16)  can  be  easily  solved  if  we  know  analytical  expressions  for  downward  and  upward  irradiances  and 

connection  between  average  cosine  [A  and  inherent  optical  properties  Cl  and  bB.  By  integrating  Eq.  (16)  over  downward 
and  upward  parts  of  a  solid  angle,  we  obtain  the  following  system  of  equations  for  downward  and  upward  irradiances: 
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Lik(z)Ek(z)  =  fi(z ),  i - 1,2, 


(17) 


where  index  1  corresponds  to  a  subscript  r/,  index  2  corresponds  to  a  subscript  u,  a  summation  is  assumed  over 


repeated  indices,  /j(z)  =  f2(z)  =  f(z)  =  2nQ(z),  and  the  differential  matrix  operator 

Lik  is  given  by  the  following 

equation 

—  4-  (9  —  n)(n  +  h  _ 

_ —(9.  4-  n)h - 

— - — T  - p'A'  1  B  * 

dz 

\  r  '  a 

d 

(18) 

-(2-Ji)bB 

-—  +  (2  +  n)(a  +  bB) 
dz 

with  the  average  cosine  / u  defined  through  the  inherent  optical  properties  a  and  bB  as  follows: 


a  — 

- ,  — ,  or  jU 

a  +  3bB  +  -y jbB  (4  a  +  9  bB) 

where  g  =  bB  /(a  +  bB  )  is  a  Gordon’s  parameter. 

In  order  to  solve  Eqs.  (17)  and,  subsequently,  Eq.  (16)  we  need  to  formulate  boundary  conditions.  Because  we 
assume  that  both  upper  (sea  surface)  and  lower  (sea  bottom)  boundaries  are  Lambertian  reflectors,  we  define  the 
conditions  as  follows: 

Eu(zb )  =  As  [Ed(zB)  +  E](ze)l  Ed( 0)  =  As  Eu( 0),  (20) 

where 

E](zb)  =fJdcpf'ond[iLq(zB,w)  =  2 nLODfo  e^^ndn  +  Los  (21) 

is  a  downward  irradiance  by  the  light  of  the  source  at  the  bottom.  The  boundary  conditions  (20)-(2 1 )  will  be  used  for 
Eqs.  (17)  as  well  as  Eq.  (17). 


1  ~g 


\l  +  2g  +  Jg(4  +  5g) 


(19) 


3.  SOLUTION  OF  RADIATIVE  TRANSFER  PROBLEM  FOR  VERTICAL  IRRADIANCES  OF 

LIGHT 

Let  us  solve  Eqs.  (17)  for  vertical  (downward  and  upward)  irradiances  of  underwater  diffuse  light  with  the 
boundary  conditions  given  by  Eqs.  (20)-(21).  The  complete  solution  to  this  problem  is  the  sum  of  general  and  partial 
solutions  to  Eqs.  (17).  The  general  solution  is  a  solution  to  Eqs.  (17)  with  the  right  side  equal  to  zero.  The  partial 
solution  is  expressed  through  a  Green’s  matrix  of  Eqs.  (17).  The  Green’s  matrix  Gik(z )  is  a  solution  to  the  following 
matrix  differential  equation: 


Lik(z)Gkl(z)  =  6(z)6il,  i  =  1,2,  (22) 

where  8tl  is  a  2x2  unity  matrix  (or  Kronecker’s  symbol).  It  is  easy  to  show  [1]  that  the  Green’s  matrix  is  expressed  as 
follows: 


1  *0 

H(z)e-a~>  | 

R0R» 

R0  H(-z)ea°z 

K  R0R co 

Q-R0RJ 

R„ 

1  (1-^co)’ 
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and  H(x)  is  a  Heavyside’s  or  step  function: 


1,  xaO, 

n  (26) 

[0,  x  <  0. 

Using  Green’s  function  (23)  we  can  write  the  solution  to  Eqs.  (17)  as  follows: 

1  R 

E,(z)=  +  Dniea°z  +  fj  Gik(z- z')fk(z')dz' ,  *  =  1,2,  m.  =  ,  ni  =  °.  (27) 

'*'‘00  1 

By  inserting  Eq.  (23)  into  Eq.  (27)  and  introducing  the  following  three  auxilary  functions, 

g(z)  =  2 Jtf0  Q(z')ea~z  dz' ,  h(z)  =  2 Jtf  Q^e^'dz' ,  m(r)  =  2jVr">  dp, ,  (28) 

we  can  simplify  Eqs.  (27)  and  obtain  the  following  equations  for  downward  and  upward  diffuse  irradiances  of  light  in 
shallow  water  body  or  arbitrary  turbidity: 


1 

,  n.  - 

R« 

*  i 

1 

E  Az)  =  A  e~a'z  +  DR0ea°z  + 


l  -  R0R« 


Ell(z)  =  ARae-a*z+Dea°z  +  Ri 


i  -R0R» 


1  +  R 

- l8(z)  +  R0  p  Z  ea°‘h(z), 

1  ^0^00 

(29) 

1  +  R 

e-a'zg{z)  +  —-^ea°zh(z). 

(30) 

The  two  coefficients  A  and  D  in  Eqs.  (29)-(30)  are  fully  determined  through  the  two  boundary  conditions  and  can  be 
written  as  follows: 

A  =  Klmh(0HDl  '  CO 

(1  -R„A,)e-"-K(A,-R„)e--‘- 

where 

As  —  R0  1  +  R0 

K  ==  — - co  = - — ,  (32) 

i  -  asr„  \-R0R«  . 

and 

E«(z)  =  nL0Dm(ccz)  +  p0L0Se-az,f,°  (33) 

is  a  downward  irradiance  by  external  sources  (direct  and  diffuse  solar  light  and  diffuse  light  of  the  sky  penetrated  into  the 
water  body). 
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4.  SOLUTION  OF  RADIATIVE  TRANSFER  PROBLEM  FOR  RADIANCES  OF  SCATTERED 

LIGHT 


Equations  (9)-(10)  for  radiances  can  be  rewritten  as: 


d 

+  « 

dz 


QAz), 


(34) 


QAz )  =  ^|^[(2  -  &EAz)  +  (2  +  mSz)  +  Ks  ^z/"°]+  lod  q(<*z)  J,  (35) 

where  Ed  (z)  and  Ed  (z)  are  given  by  Eqs.  (29)-(33). 

With  known  right  side  the  solution  of  Eq.  (34)  for  scattered  light  radiances  can  be  obtained  as  a  sum  of  partial 
and  general  solutions: 


L(z,fi)=f0Bdz'QAz')G(z-z',v)  + 


f  C,eaz!\ 
\C2e-azl\ 


iU&O, 

A*  <  0, 


(37) 


where  G(z,fl)  is  a  Green’s  function  of  Eq.  (34)  that  satisfies  the  following  equation: 


+  a]  G(Z,A)  =  6(z), 
dz  1 


(38) 


and  is  expressed  as  follows: 


G(z,n)  =  H 


z)e 


■a  z/fi 


\l)  |^| 


(39) 


where  HQ c)  is  a  Heavyside’s  function  defined  by  Eq.  (26).  The  insertion  of  Eq.  (39)  into  Eq.  (37)  gives  the  following 
result. 


L{z,y)  = 


e-°z/ft  [c,  +fo  QAz’)eaz'^  dz'/fij,  s  0, 

[c2  +  £  QAz')e-a  Z'M  dz'/\  H  |],  A*  <  0. 


The  constants  C,  and  C,  are  determined  by  the  boundary  conditions  given  by  Eqs.  (20)-(21), 


c,  —  (C2  +  £q)>  c2  Ab 


£n  +  As  £o  +  Lon  +  Vo  (W”)g~°  ^ 

Asp) 


where 


s0  =  2£"  Qt (z) n(a z) dz,  eB  =  2£"  QT (z) n[a( zB  -  2)] dz. 


(40) 


(41) 


(42) 
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(43) 


m_  =  m(azB )  =  e'azn^  >  0,  m+  =  m{-azB )  ■  2 J*ojud/4  e“v''  >  0,  n(z)  =JV"^  dpi .  (44) 

Equations  (40)-(44)  fully  define  scattered  light  radiances.  By  adding  direct  radiance  (7)  we  obtain  the  total  radiances  in 
the  shallow  sea:  _ _ _ 


L,(z,H,<p)  =  \ 


e-az/K 


+  C,  +  foQT(z')eaz‘/t‘  dz'fri ]  +  Lose-azltl"  d(fi-  pi0)S((p- it), 
e~°Zh‘  [C2  +fz“  &(*')*“" lVW  dz'/\  pi  |  , 


jU  &  0, 

(45) 

pi<o. 


Equation  (45)  represents  a  self-consistent  solution  to  the  radiative  transfer  problem  for  radiance  of  light  in  shallow  water 
with  wavy  Lambertian  surface  and  Lambertially  reflecting  bottom.  It  takes  into  account  multiple  scattering  inside  water 
column  and  from  the  bottom  and  the  surface,  and  it  is  valid  in  the  whole  range  of  variability  of  inherent  optical 
properties  (0  £  <2  £  °°,  0  <  /}  £  00 ,  0  £  Z)  <  0.5)  and,  consequently,  can  be  used  for  waters  with  arbitrary  turbidity 
and  absorption.  For  modeling  purposes  inherent  optical  properties  can  be  generated,  for  example,  using  approach 
proposed  in  Ref.  [5]. 


5.  HORIZONTAL  IRRADIANCES  OF  LIGHT  IN  SHALLOW  WATERS 

The  horizontal  irradiance  in  shallow  water  is  a  function  of  depth  and  azimuth  angle  (p  between  the  direction  of 
solar  rays  and  the  direction  of  light  we  measure.  It  is  given  by  the  following  equation: 

Eh (z,(p)  =  n |[Lod  +  C,]<E>(a z)  +  C2 0(-a z)  +J*“  W(a[z- z'])Qr(z')dz'  J 


here 


J Los  ea  zlf,°  cos (*:  -  <p),  |  x  -  (p  |<  n  /  2, 

j  0,  \n  -q)\>n  12. 

W(r)  =  -ite***  dpi/ pi. 


(46) 


(47) 


the  function  W(t)  has  an  infinite,  narrow  and  integrable  peak  at  T  =  0,  so  it  can  be  replaced  by  a  delta  function: 


W(t)s^6(t),  <5(az)  =  —  d(z). 
4  a 


(48) 


This  significantly  simplifies  the  expression  for  horizontal  irradiance: 

Eh  (z,cp)  =  x  |[Lod  +  C,  ]<&(a  z)  +  C2  $(-a  z)  +  QT  (z)  J 

+  \L0S^\-  pil  zn‘°  cos (ji  -cp),  \n-cp\<.n  1 2, 
|  0,  \x  -  cp\>n  12, 


(49) 
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with  the  coefficients  C,and  C2  given  by  Eqs.  (41).  It  is  necessary  to  note  that  the  third  component  in  the  right  part  of 
Eq.  (49)  is  not  equal  to  zero  only  on  a  sunny  side  of  the  viewing  direction  \k  -  (p  \&7t  12. 


6.  CONCLUSION 

We  obtained  here  a  simple  self-consistent  solution,  given  by  Eq.  (45),  to  the  radiative  transfer  problem  in  a 
shallow  sea  with  wavy  surface  and  diffusely  reflecting  bottom.  This  solution  takes  into  account  multiple  scattering-of 
light  inside  water  body  as  well  as  multiple  reflections  from  bottom  and  wavy  surface.  It  is  valid  for  waters  with  arbitrary 
absorption  and  turbidity.  The  obtained  solution  was  used  to  calculate  horizontal  azimuthally  dependent  irradiances  in  the 
shallow  water  with  reflecting  bottom  and  wavy  surface.  Expression  (49)  for  horizontal  light  irradiances  is  a  principal 
equation  that  will  allow  to  solve  the  problem  of  horizontal  (or  diver)  visibility  of  an  object  submerged  in  a  shallow 
water  with  arbitrary  inherent  optical  properties. 
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